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We propose a new method to test the (resummed) next-to-leading- 
order BFKL evolution kernels using the Mellin transformed j-moments of 



, the proton structure function F2 ■ 

^ ' 1. How to test BFKL evolution equations? 

£^ [ The Balitsky Fadin Kuraev Lipatov (BFKL) evolution equation [Q, de- 

^5 rived in the framework of perturbative QCD, has held the attention of the 

scientific community since a long time. The summation of leading loga- 
rithms of energy in the perturbative expansion gives valuable tools for the 
investigation of deep-inelastic scattering at small XBj (equivalently large 
c I ! energy squared W 2 ~ Q 2 /xsj)- Indeed, the first experimental results from 

HERA confirmed the existence of a strong rise of the proton structure func- 
tion F2 with energy in agreement with the trends implied by the solution 
r> I of the BFKL equation. It has been possible [2] to describe the old data at 

small XBj and in a certain range of Q 2 . However the price to pay was to 
get a phenomenological value of the intercept (the exponent of 1/x Bj in the 
BFKL formula, see later in the text) less than the predicted range (the cor- 
responding value of the strong coupling constant is a s ~ .1 instead of ~ .2). 
This was revealing the need for rather sizable higher order corrections. 

At the next-to-leading log level, these corrections have been calculated 
after much efforts [3] and appeared to be so large that they overshoot the ex- 
pected phenomelogical effect and could even invalidate the whole approach. 
Soon after, it was realized [H that the main problem came from the exis- 
tence of spurious singularities which ought to be cancelled by an appropriate 
resummation at all orders of the perturbative expansion. This is required 
by the QCD renormalization group. Indeed, various resummation schemes 
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have been proposed [U El E] which satisfy the renormalization group re- 
quirements while retaining the exact value of the next-leading log term in 
the BFKL kernel computed in Refs.0. Hence, the constraints can be satis- 
fied and the next-to-leading order introduced without destroying the whole 
scheme. 

However, the resummation schemes possess some ambiguity, since higher 
order logs (beyond the next-to-leading ones) are not known. Such variations 
appear, e.g. in Ref.|3], where four different resummation schemes have been 
proposed. These schemes, denoted scheme l,..-4 m the following, will be 
the subject of our present study. Other schemes have been proposed El 
and will be studied as well later [7j. It is worth to confront these various 
schemes with data, to check their validity and distinguish between different 
resummation options. 

Precise phenomenological tests of QCD evolution equations are one of 
the main goals of deep inelastic scattering phenomenology. For DGLAP 
evolution 8 , it has been possible to test it in various ways with next-to- 
leading log Q 2 (NLO) corrections and it works quite well in a large range of 
Q 2 . Testing precisely BFKL evolution beyond leading order is much more 
difficult. The main problem is the complicated mismatch between QCD per- 
turbative and non perturbative inputs, since the corresponding factorization 
properties, i.e. kx factorization [5], are more involved than for DGLAP evo- 
lution. A way out could be to stay within the perturbative regime by using 
only massive or highly virtual colliding particles, like 7*— 7* scattering, but 
the data are yet too imprecise that no definite conclusion can be drawn. 
Note also that some perturbative QCD ingredients (such as the so-called 
"impact factors") are not yet but will be soon available |lUj . 

In the present paper, we propose a method for testing the (resummed) 
BFKL predictions for the proton structure functions, via a transformation 
to Mellin space. 

On the one hand, the present set of data allows for a precise determina- 
tion of the Mellin transform of F2 in a large range of Q 2 and j, the Mellin 
conjugate of xbj, considered as a continuous variable. On the other hand, 
the BFKL predictions at leading order and beyond are easier to formulate 
in Mellin space since one can obtain tests of the evolution kernels which 
are essentially dependent on the calculable perturbative part. In our for- 
mulation, we assume that kx factorization remains valid for the resummed 
kernel. An improved formulation of kx factorization containing the NLO 7* 
impact factors when they will be available, will allow to refine our study in 
the future. 

The proposed method has the following features. It treats in paral- 
lel both LO and (resummed) NLO BFKL kernels. It uses a kx factorized 
formulation of the structure functions which includes the factorized Green 
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function solution derived in Ref. (Sj. The present essay is of introductory 
nature, and present a first, non sophisticated, phenomenological investi- 
gation of the proton structure functions [111 I12j . where we compare the 
LO-DGLAP GRV parametrization |13j of proton structure functions with 
the BFKL kernel predictions. In a forthcoming publication [7], we shall give 
an extensive and more systematic study using our method. 

The plan of contents is organized as follows; In the next section, we 
start by expressing the LO BFKL predictions in Mellin space, defining three 
characteristic relations. In section 3 we elaborate the corresponding set of 
predictions for the resummed NLO-BFKL kernels. An application to the 
DGLAP/BFKL comparison is presented in section 4. Section 5 provides a 
conclusion and an outlook on future work. 



2. BFKL predictions in Mellin space 

The formulation of the proton structure functions in the (LO) BFKL 
approximation can be expressed as follows |2j: 

where one has written the BFKL kernel as 

Xlo(7) = 2V(1) - V>(7) - ^(1 - 7) • (2) 

In formula with conventional notations, Ft,Fl,G stand respectively 
for transverse , longitudinal and gluon structure functions, a s is the (fixed) 
coupling constant, a; (7) is an (unknown) non-perturbative coupling to the 
proton while 

(h T \_ a s (r(l-7)r(l + 7 )) 3 1 /(i +7 )(i_2)\ 

\h L J 3tt 7 r(2 - 2 7 )r(2 + 2 7 ) 1 - | 7 ^ 7(1-7) )' { ) 

correspond to the known perturbative couplings to the photon, usually 
called LO "impact factors" in the literature. Note that, in the framework 
of kx factorization, 7 plays the role of a "running" anomalous dimension, 
whose physical value is determined by the integration of (^Q). As already 
mentionned in the introduction, formula (^Q) gives rise to an interesting ef- 
fective BFKL phenomenology [2] in the small x Bj region, but it has to rely 
on the parametrization of the unknown non perturbative function ^(7) and 
leads to values of a s quite smaller than expected. 
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Mellin-transforming in j-space, one easily finds 



Fl I = / \ ^ I ■ ; n . N . — I h L | w( 7 ) • (4) 

G 



Looking for the poles in 7, it is straightforward to use the residue formula 1 
and get 

where 7i(j) are the (7$ < 1/2) roots of the equation 

J - 1 = ^-^XLohiij)) ■ (6) 

7T 

In fact, with good accuracy at large enough Q 2 , (comparable to the leading 
twist approximation in DGLAP evolution), one can only retain the right- 
most pole 71 (j). We are thus left with the following simple formula as a 
starting point of our analysis: 

Ft\ 1 /O 2 \ 7l(j) /M7it?'))\ 

l M T 0)) J u <71(J,) - (7) 

From equation (J7J), three model- independent predictions, i.e. indepen- 
dent of non-perturbative assumptions, can be drawn: 

i) The Mellin transform of Fz = Ft + Fl should verify: 

ln J F 2 (j,Q 2 ) = 7l (j)ln(Q 2 ) + /(j) (8) 

in some range of j near 1 where the BFKL equation is expected to be 
relevant. The function f(J) regroups all Q 2 -independent terms in l[7)). 

ii) 71 (j), extracted from (jJJ) as the slope in ln(Q 2 ), should verify the 
equation @ for the anomalous dimension, namely 

XLodiU)) = -V (j-i) , (9) 

with constant ct s , and xlo given by ©. 



Care is to be taken of the contour at infinity, see the second reference of |7|. 
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iii) The gluon structure function (one may also choose the obervable 
Fl) should verify, via Mellin transform: 

ln(GU Q 2 )) = In (F 2 (j, Q 2 )) - 7l (j) In (h T + h L ) . (10) 

The predictions (jSJ), © (HOj) represents a stringent set of constraints which 
have to be verified by the Meflin-transformed of the proton structure func- 
tions in a region j near 1. In fact we will confirm that the (LO) BFKL 
kernel does not pass this step. 



3. Resummed NLO-BFKL predictions in Mellin space 

Interestingly, using a reasonable kx factorized ansatz 2 , the predictions 
i)— iii) remain valid for NLO-BFKL resummation kernels, up to specific mod- 
ifications due to the running of the coupling constant. 

Let us formulate the (resummed) NLO-BFKL structure functions in 
Mellin space as follows: 



^ X ^ h L rj( 7 ) , (11) 





where, by definition 

— X( 7 ,j)= X NL0(l,j)- (12) 

The function X{^i,j) appears in the solution of the Green function de- 
rived 3 from the renormalization-group improved small-x^j equation jS], 
XNLoild) is the resummed NLO-BFKL kernel and b = 11 - 2/3 N f /N c 
defines the running of the coupling constant 

- *s(Q 2 ) = s • (13) 

bin (Q /A-qcd) 

Before going further, let us comment formula (|lljl - This formula cap- 
tures the (large) Q 2 -dependent part of the gluon Green function which has 
been shown to have factorization properties 5 . In fact the non perturbative 
contribution has been factorized out in the function 77(7). Some unknown 



2 kr factorization has not been yet proven at NLO, but is a reasonable ansatz fulfilling 
the known theoretical requirements on the kernel properties discussed in 0. 

3 The second variable of X(<y,j) in 1 1 2t corresponds to the choice of a reference scale 
fj, — > j — 1 dictated by the treatment of the Green function fluctuations near the 
saddle-point 5 . 
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(^-dependence may still remain in the NLO contributions to the impact 
factor vector (7t), which we neglected in the present analysis. 

Starting from this ansatz, let us derive the NLO constraints similar to 
(|5])-(|1U[). At large enough Q 2 /Aq CD , one can use the saddle-point appox- 
imation to evaluate Assuming that the perturbative impact factors 

and the non-perturbative function r] do not vary much 4 , the saddle-point 
condition reads 

• i 1 N c a s (Q 2 ) 

J-l~ TZTT~k V XNLO{l,3) = XNLO{j,3) , (14) 



b ln(Q 2 /A 



QCD 

where j(j,Q 2 ) is the saddle-point value. Relation (|14|) is nothing else than 
the NLO extension of condition ii) of Q to the case of a running coupling 
constant JT3J). 

Inserting the saddle point defined by (|14j) in formula one obtains 
the new set of constraints at NLO level as follows: 
i) The Mellin transform of F<i = Ft + Fl verifies: 

■lnF 2 (j,Q 2 ) ~ 7(J,Q 2 ) (15) 



d\n(Q 2 



where Q 2 ) is now a smoothly (^-dependent effective anomalous dimen- 
sion defined by the following property. 

ii) 7(j, Q 2 ) verifies the anomalous dimension equation, namely 

XNLo(l(j,Q 2 )) = j^r T (j-1) , (16) 

where xnlo is ° ne of the resummed NLO-BFKL candidate kernels proposed 
in the literature. 

iii) The gluon structure function (one may also choose the obervable 
Fl) verifies, via Mellin transform: 

ln(G(j,Q 2 ))=ln(# 2 (i,Q 2 ))- 7 (j,g 2 ) In [h T ft) + h L {l)} , (17) 

where NLO effects of impact factors are neglected. 

The interest of the relations Q15|) H16|) ()17|) is that they are formally 
similar with the LO ones by the direct substitution of the LO kernel by the 
NLO ones and of a fixed coupling constant by the running one at one-loop. 



4 We do not take into account modifications e.g. coming from powers of 7 in the 
prefactors which may shift the saddle point We thus assume a smoothness of 
the structure function integrand around the saddle-point in agreement with the phe- 
nomenology Q. 
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They are only approximate, since they rely on a saddle-point approximation 
which may not be always justified (see for a discussion). However, in 
the present context, the validity of the saddle-point approximation can be 
tested directly from the phenomenological analysis. Due to the observed 
smoothness of the Mellin transforms, we do not expect large corrections to 
the saddle-point results. 

4. Application: the "proximity" between DGLAP and 

NLO-BFKL 

In this section, we want to check the reliability of the Mellin space 
method by a study of DGLAP parametrizations of the data . It is 
well-known that DGLAP parametrizations fit well the data in a large range 
of XBj and Q 2 > 1 GeV. Choosing such a parametrization of structure 
functions, namely the GRV set of structure functions JB], we are able to 
Mellin transform them easily, and thus discuss the comparison between 
DGLAP and LO/NLO BFKL evolution equations. DGLAP evolution is 
automatically obeyed by the input functions and we want to compare them 
with BFKL evolution using relations © © dB for LO and {EJ CEJ) (JUJ) 
for NLO. The physical question we ask in this application is whether or 
not there may exist a compatibility between DGLAP and BFKL evolution 
equations. 

Let us first consider the singlet density distribution in Mellin space £ = 
(q+q)(j), see Fig.l. It is an easy exercise to obtain it from the input GRV 
parametrizations and LO DGLAP matrix elements in Mellin space (Hj. 

It is clear from Fig.l that there exists an interval 1.3 < j < 1.7 in which 
the slope of In X as a function of InQ 2 is almost constant. The observed 
approximate constancy meets the requirement 5 of condition i). We will 
focus our study to this region which is Mellin-conjugated to the small xsj 
region. 

Taking into account the anomalous dimension values 6 , it is now straight- 
forward to look for the LO prediction @ and the various NLO predictions 
1|16|) depending on the choice of resummation scheme in [I]. 

We have displayed in Fig. 2 for comparison the results for both the 
standard (LO) BFKL and for one of the Resummed NLO BFKL schemes 
{scheme 4)- In our application plotting x(7i(j«)) as a function of ji should 
give points aligned on a straight line extrapolating to at j = 1. The slope 

5 The NLO condition l|15|l implies a smooth variation of the slope. Presently, we will 
study the Q 2 average only, delaying a more refined (but deserved) study of the Q 2 
dependence of the slope for a forthcoming publication jjj. 

6 In this preliminary study we considered only discrete values Ji(ji), where ji = 
1.3, 1.4, 1.5, 1.6, 1.7. determined for a fixed range 2 < log Q 2 < 6. 
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Fig. 1. Evolution of the Mellin transformed (q + q)(j) as a function of Q 2 . 
From left to right and top to bottom: j = 1.1, 1.3, 1.4, 1.5, 1.6, 1.8. 



gives the (average) value of nNc/a s . By simple inspection of Fig. 2, it is 
clear that there is a large difference between LO and NLO (scheme 4 of 
0]) results. The LO test completely fails in shape and magnitude, while 
the NLO test is satisfactory. In this case, the measured slope leads to an 
average value a ~ .23 which is reasonable for the Q 2 range considered for 
the structure functions. 

Some comments are in order. We have here used a set of structure 
functions which, on the one hand verifies the DGLAP evolution and, in 
the other hand give a satisfactory fit of data. The conclusion of our test 
would be that there is not much difference between DGLAP evolved effective 
anomalous dimensions from DGLAP evolution and from some well-choosen 
resummed NLO-BFKL kernels. This remark corroborates the proximity of 
DGLAP and NLO-BFKL predictions for the gluon anomalous dimension in 
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Fig. 2. Study of the LO and NLO BFKL kernels in Mellin j-space. 
Top: (LO) BFKL kernel; 

Bottom: Rcsummed NLO BFKL kernel, scheme 4 of Rcf. 0]. 
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[2j and, in a different context, the smallness of BFKL-like corrections to the 
DGLAP evolution equation found in |14j for HERA data. 

The new point is that our method allows for sensitive tests of the re- 
summed NLO-BFKL kernels. For instance the schemes 1,2 fail ^21 an d this 
is to be related to the fact that not all large logs are resummed @J. Scheme 
3 gives a satisfactory shape but with an averaged a s ~ .27 quite too strong. 
This is for illustration of the sensitiveness of the method. 

We postpone to a more systematic study jZj the comparison between 
different sets of structure functions, different methods, and also the consid- 
eration of non-DGLAP parametrizations in order to evaluate the systematic 
biases which could occur from these different options. We also leave for fu- 
ture study the relations (|10 j l.p7j ) which gave promising results in |llj . 

5. Conclusion and outlook 

• We have proposed a method for confronting with precision "data" 
the various resummed BFKL kernels with next-to-leading log accu- 
racy. These "data" are the Mellin-transformed of the proton structure 
function. 

• Due to the high precision of modern experimental data on F%, we 
expect the Mellin transform to be well determined, at least in the 
region of j and Q 2 needed for the test. 

• We make use of a kx factorized formulation of the structure function 
which grasps the constraints coming from the QCD renormalization 
group improved small-xsj equation [3]. 

• In a first, preliminary, application using the parametrization ^3] as 
input, the method leads with high sensibility to an incompatibility 
of the LO BFKL kernel while only one of the four resummed NLO 
BFKL schemes of Ref. [3J shows a compatibility with the DGLAP 
parametrization of data. 

The application we have performed is far from complete, and was only 
serving as a test of the method sensitivity. Further studies in various inter- 
esting directions deserve to be pursued . let us quote some of them. 

It is interesting to test whether the present set of data and with which 
accuracy, the Mellin transformed of the structure functions can be obtained. 
In particular, it would be useful to see eventual differences between DGLAP 
and non-DGLAP parametrizations of data, to see whether and where there 
could be a systematical bias introduced by the DGLAP framework. 
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Concerning the resummation schemes, the application of the method 
with the correct NLO accuracy requires to take fully into account the cou- 
pling constant running. It is thus required to separate data in small regions 
of Q 2 and make the tests separately in each region, in order to follow the 
Q 2 evolution of the coupling constant. 

We expect to be able to answer these questions soon (7j. 



The materials for the phenomenological application comes from Julien 
Salomez jll|H2j . whom I thank. I thank Gavin Salam for fruitful discussions 
and Christophe Royon and Laurent Schoeffel for a fruitful collaboration. 

As well as many of my colleagues, I am grateful to Jan (Kwiecinski) for 
his enthousiasm and readiness of mind, and for his work which represents a 
cornerstone in our domain. Bonne et Longue Continuation, Dear Jan! 



[1] L.N.Lipatov, Sov. J. Nucl. Phys. 23 (1976) 642; V.S.Fadin, E.A.Kuraev 
and L.N.Lipatov, Phys. lett. B60 (1975) 50; E.A.Kuraev, L.N.Lipatov and 
V.S.Fadin, Sov. Phys. JETP 44 (1976) 45, 45 (1977) 199; I.I.Balitsky and 
L.N.Lipatov, Sov.J.Nucl.Phys. 28 (1978) 822. 

[2] H Navelet, R.Peschanski, Ch. Royon, S.Wallon, Phys. Lett. B385 (1996) 357. 
S.Munier, R.Peschanski, Nucl.Phys. B524 (1998) 377. 

[3] V.S. Fadin and L.N. Lipatov, Phys. Lett. B429 (1998) 127; M.Ciafaloni, Phys. 
Lett. B429 (1998) 363; M. Ciafaloni and G. Camici, Phys. Lett. B430 (1998) 
349. 

[4] G.P. Salam, JEEP 9807 (1998) 019 

[5] M. Ciafaloni, D. Colferai, G.P. Salam, Phys. Rev. D60 114036, , JEEP 9910 
(1999) 017; 

M. Ciafaloni, D. Colferai, G.P. Salam,A.M. Stasto, Phys. Lett. B541 (2002) 
314. 

[6] Stanley J. Brodsky, Victor S. Fadin, Victor T. Kim, Lev N. Lipatov, Grigorii 
B. Pivovarov, JETP Lett. 70 (1999) 155. 

[7] R.Peschanski, Ch. Royon, L.Schoffel, to appear. 

[8] G.Altarelli and G.Parisi, Nucl. Phys. B126 18C (1977) 298. V.N.Gribov and 
L.N.Lipatov, Sov. Journ. Nucl. Phys. (1972) 438 and 675. Yu.L.Dokshitzer, 
Sov. Phys. JETP. 46 (1977) 641. For a review, see e.g. Yu.L. Dokshitzer, 
V.A. Khoze, A.H. Mueller, S.I. Troyan, Basics of Perturbative QCD . 

[9] S. Catani, M. Ciafaloni, F. Hautmann, Nucl.Phys. B366 (1991) 135; S. 
Collins, R.K. Ellis, Nucl.Phys. B360 (1991) 3; E.M. Levin, M.G. Ryskin, 
Yu.M. Shabelski, A.G. Shuvaev, Sov.J.Nucl.Phys. B53 (1991) 657. 



Acknowledgement 



REFERENCES 



12 



BEYOND 1 PRINTED ON FEBRUARY 1, 2008 



[10] J. J. Bartels, D. Colferai, S. Gieseke, A. Kyrieleis, Phys.Rev. D66 094017, and 
references therein. 

[11] J. Salomez, Le modele des dipoles en QCD perturbative (in French), Saclay 
preprint T02/147 (2002), Di ploma Memoir for the "P EA Rhne-Alpin, ENS 
Lyon, France". Available at: http://www-spht.cea.fr/articles/t02/147/ 

[12] J. Salomez, unpublished notes. 



[13] M. Gluck, E. Reya, A. Vogt, Z.Phys. C67 (1995) 433; See M. Gluck, E. Reya, 
A. Vogt, Eur.Phys.J. C5 (1998) 461, for updated parametrizations. 

[14] R.S. Thorne, Phys.Rev. D60 (1999) 054031;G. Altarelli, R.D. Ball, S. Forte, 
Nucl.Phys B621 (2002) 359, and references from the same authors therein. 



